Thermalization of a particle by dissipative collisions 
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One considers the motion of a test particle in an homogeneous 
fluid in equilibrium at temperature T, undergoing dissipative colli- 
sions with the fluid particles. It is shown that the corresponding lin- 
ear Boltzmann equation still posseses a stationary Maxwellian veloc- 
ity distribution, with an effective temperature smaller than T. This 
effective temperature is explicitly given in terms of the restitution 
parameter and the masses. 

PACS numbers: 51.10.+y, 44.90. +c 



The search for stationary states of granular matter has recently been a subject 
of interest for both experimental and theoretical reasons [1,2]. Granular matter can 
be modelized by spherical particles that partially dissipate their kinetic energy at 
collisions. If (u, v) and (u, v) denote the velocities of two spheres of mass m and 
M before and after the collision, they are related by 

mu + Mv (1) 
-cmt-(v-u), 0<a<l (2) 
ct X • (v - u) 

where <x is a unit vector normal to the surface of the spheres at the point of 
impact, and cr 1 ' points in the tangent direction: & ■ & = 0. The first relation is 
the conservation of the center of mass momentum, whereas the second one says 
that the normal component of the relative velocity reverses its direction with a 
magnitude reduced by the factor a, the so called restitution parameter (the tangent 
component remains unchanged). Solving (|l]), (^) for u and v gives 

u = u+(l-/i)(l + a)(d--(v-u))6-, v = v - fi(l + a)(a • (v - u))<7 (3) 

with fj, = m/(m + M). The inverse relation is obtained by exchanging the roles of 
initial and final velocities and a into oT 1 in (|3|). 

We consider now the motion of a single particle of mass M (the test particle) 
in a fluid of particles of mass m in thermal equilibrium at temperature T by means 



mu + Mv = 

<T ■ (v — u) = 

a 1 - ■ (v - u) = 



1 



of the linearized Boltzmann equation. All particles are spheres of diameter d and 
the test particle undergoes inelastic collisions characterized by the restitution pa- 
rameter a with the particles of the host fluid. Then, in the homogeneous situation, 
the probability density /(v, t) for finding the test particle at time t with velocity 
v obeys the kinetic equation 



0_ 

Of 



/(v, t) = pd 2 J duj d&(& • (v - u))0(<7 • (v - u)) x 



[a- 2 f (v - n (1 + a" 1 ) (<r • (v - u))<x, t) <fr (u + (1 - ft) (l + cT 1 ) (a • (v - u))<x) 

-/(v,t)0r(u)] (4) 

Here p<pr(u) is the equilibrium state the fluid having uniform density p and 
Mawellian velocity distribution 

pm \ I mu \ 1 



Mu)=^j exp^/3— j, „= H) /3=— (5) 

In (^|), u is the velocity of a particle in the fluid, v is the velocity of the test 
particle and the <r-integral runs on the unit sphere. The velocity arguments of / 
and 4>t in the gain part of the collision term in (Q) are precisely those that the 
particles must have before the collision to produce post-collisional velocities u and 
v, according to (||). The term(<r • (v — u))0(<r • (v — u)) (with 9(y) the Haeviside 
function) reflects the fact that the collision frequency is proportional to the velocity 
of approach of the colliding pair (free motion between encounters). The factor or 2 
compensates for the contractive transformation (^|) in velocity space and ensures 
that the normalization of / is properly conserved in the course of the time. 

The point of this letter is to show that, despite of its dissipative collisions with 
the particles of the bath, the distribution of test particle still posseses a stationary 
Maxwellian velocity distribution, with a reduced effective temperature T that 
will be explicitly given in terms of the restitution parameter and the masses. 

We look for a stationary state setting simply /(v, t) = /(v) and equating the 
r.h.s. of mj to zero. With the change of variable w = v u, / has to satisfy 



= J da- J dw(a ■ w)9{a ■ w) x 



[a 2 f (v - p (I + a x ) (a ■ w)<r) (f) T (v - W + (1 - p) (l + a x ) (a ■ w)<r)) 

-/(v)0r(v-w)] (6) 
Let us perform the w— integration at fixed <r in a frame 

w = wi&i + w 2 a 2 + w 3 a-3 (7) 
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where {ar± = cr, a%, ^3} is an orthonormal system of unit vectors with &s in 
the plane perpendicular to er. Note that the argument of / does not depend on the 
variables wi and W3, so the corresponding integrals can be carried out directely 
on 4>t- In terms of the variables (Q) we have 

(j>T (v — w + (1 — /i) (l + a -1 ) (o - • w)<r)) 

f3m\ 3 ^ 2 ( f3m ,„ , , , _ n , ,2] 

ItT j CXP 1 2 _ ( <T ' v+ ( a -M(! + a jj^ij f 

/ (3m. „ , 2 \ / /3m. „ , 2 \ . . 

x exp — — (<r 2 -v-tua) exp — — (<r 3 • v - to 3 ) (8) 



Hence the integrations on 1U2 and u>3 can be readily performed here (as well in 
0t(v — w)) and the stationary equation takes eventually the form 



0= da dyy 



/(v ;/'/<r) - xi) I — ■ v + (1 - ?i)y) 2 



-/(v.) exp I • v - y) 2 



(9) 



To obtain (Q) we also put »i = ay in the first part of the integrand in (g) (gain 
term) and we defined the new parameter 7/ = fi(a + 1). 
Now, let us set 

/(v) = exp (-/?Af7j 2 /2 7 ) , « = |v| 

omitting the normalization factor. One notices then that the following equality 
holds 

Af7 -1 |v - r]ya-\ 2 + m(& • v + (1 - rj)y) 2 = Mj~ 1 v 2 + m{a ■ v - y) 2 (10) 
provided that 7 -1 has the value (independent of v, <r and y) 

1= m(2-r,) 
' Mr) K ' 

The relation ( |l0|) implies vanishing of the integrand in (Q) . Introducing the effective 
temperature by 

T cS = jT (12) 

we conclude that the Maxwellian /(v) = <fi T ca(y) is a stationary distribution for 
our test particle. The explicit formula for the factor 7 reducing the temperature 
in terms of the restitution parameter and the masses reads 

a + 1 

7= 2+(l-a)m/M- 1 (13) 
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Obviously 7 = 1 when a = 1, as it should. For a fixed mass ratio m/M, 7 
decreases as a is varied from 1 to 0. Thus it takes its smallest possible value 
7 = (2 + m/M)^ 1 < 1/2 at a = when the collisions are fully dissipative. 

This result deserves several comments. Usually the dependence of the station- 
ary state on the kinematical variables is determined by the conservation laws. The 
stationary Maxwellian velocity distribution for the hard sphere Boltzmann equa- 
tion follows from the conservation of the total kinetic energy ^mu 2 + ^Mv 2 at 
collisions. It is therefore quite remarquable (and not obvious beforehand) that the 
distribution remains Gaussian when collisions dissipate the kinetic energy. Notice 
that the integrand in (^|) vanishes pointwise, that is, (f> TB n (v) is also stationnary 
in situations where the collision frequency, proportional to («r ■ ~w)9(fr ■ w), is re- 
placed by a more general positive function of <x • w. A particular case is the so 
called Maxwell gas for which the collision frequency is independent of the veloc- 
ity. Moreover, the result holds in any space dimension d = 1,2,3 with the same 
formulae (|l2|), ( |l3| ) for the effective temperature. 

These findings must be contrasted to other physical situations where the 
stationary distribution is definitely not Maxwellian. One example is the Lorentz 
model with dissipative collisions: a particle, acted upon an uniform external field, 
undergoes inelastic collisions with randomly distributed static (infinitely heavy) 
scatterers. Here the dissipation (however weak it may be) sufficies to balance the 
energy flow from the external field and to guarantee the existence of a stationary 
state. In the weakly dissipative regime, the stationary velocity distribution is found 
to behave as ~ exp(— (1 — a)v 4 /2) [3,4]. The physics of the Lorentz gas is not the 
same since there is no thermalization mechanism available for the moving particle. 
In [5] , one considers a one-dimensional system of infinitely many inelastic particles 
in contact with a thermal bath, modelized by a Fokker-Plank equation in the 
limit of weak dissipation. In this case, the stationary velocity distribution behaves 
as ~ exp(— Cv 3 ), but here dissipation occurs internally within the many particle 
system, and not in relation with its interaction with the thermal bath. 

Further investigations of interest are the study of the approach to the steady 
state in the course of time and the possible generalization of our result to a system 
of several interacting test particles. 
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